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We theoretically discuss several aspects of sequential superradiant scattering from atomic Bose- 
Einstein condensates. Our treatment is based on the semiclassical description of the process in 
terms of the Maxwell-Schrodinger equations for the coupled matter-wave and optical fields. First, 
we investigate sequential scattering in the weak-pulse regime and work out the essential mechanisms 
responsible for bringing about the characteristic fan-shaped side-mode distribution patterns. Second, 
we discuss the transition between the Kapitza-Dirac and Bragg regimes of sequential scattering in 
the strong-pulse regime. Finally, we consider the situation where superradiance is initiated by 
coherently populating an atomic side mode through Bragg diffraction, as in studies of matter-wave 
amplification, and describe the effect on the sequential scattering process. 
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I. INTRODUCTION 

The recent observation of superradiant scattering from 
atomic Bose-Einstein condensates (BECs) Q has 
opened up a new perspective on collective processes in- 
volving matter and light. In a typical superradiance ex- 
periment, an elongated BEC is exposed to a far detuned 
laser pulse with wave vector kie x traveling perpendicular 
to the condensate axis. This "pump pulse" causes some 
of the atoms to undergo spontaneous Raylcigh scattering 
and acquire recoil momentum. Due to interference, the 
recoiling atoms and the BEC atoms at rest start to form 
matter- wave gratings whose amplitudes grow rapidly in a 
self- amplifying process induced by subsequent scattering 
events. Initially, the scattered photons can go in any di- 
rection; however, due to phase matching effects induced 
by the elongated condensate shape, the fastest growth 
is experienced by those gratings for which the scattered 
photons travel up and down the long condensate axis 
with wave vectors ±ke z in the two so-called optical end- 
fire modes. Following a mode-competition process 0, 
the ensuing superradiant light emission is concentrated 
into the two endfire modes. Consequently, the recoil- 
ing atoms possess very well-defined momenta kie x ± ke z 
making 45° angles with the direction of the applied laser 
pulse, and they thus form two first-order atomic "side 
modes." To observe the side modes, one lets the system 
evolve freely for some time after the laser pulse has been 
turned off. During this time the recoiling atoms separate 
in space from the condensate at rest and become clearly 
visible as distinct spots in absorption images. 

The above scenario describes the basic mechanism of 
BEC superradiance, but the full range of effects turns out 
to be much richer. On the one hand, BEC atoms can also 
absorb a photon from an endfire mode and re-emit it into 
the laser mode. The accompanying momentum transfer 
—kie x ±ke z causes the atoms to recoil backwards at a 45° 
angle against the direction of the applied laser field. How- 
ever, this process violates energy conservation and only 
occurs for intense laser pulses shorter than the recoil time 
(Kapitza-Dirac regime). On the other hand, upon an in- 



crease in the pump-pulse duration the atoms in the side 
modes may undergo additional scattering cycles thereby 
populating second- and higher-order side modes (sequen- 
tial scattering). In the Kapitza-Dirac regime, one typi- 
cally observes an X-shaped pattern with the initial BEC 
in the center and the recoiling atoms moving both into 
and against the direction of the applied laser pulse 0. 
For weaker pulses with a duration long compared to the 
recoil time (Bragg regime), however, the side modes are 
observed to form a characteristic fan pattern involving 
forward recoiling atoms only p|. 

In recent work, we have derived a theoretical descrip- 
tion that is able to reproduce and explain these and 
other key observations of BEC superradiance in a uni- 
fied way 0, 0. Our approach which is based on the 
time-dependent Maxwell-Schrodinger equations follows 
two main concepts. First, we describe the field dynam- 
ics semiclassically. This treatment can be expected to 
be very accurate as soon as the populations of side and 
endfire modes become macroscopic and quantum fluctu- 
ations are no longer essential. The second main feature is 
the inclusion of spatial propagation effects of the matter- 
wave and optical fields along the condensate axis. It turns 
out that the consideration of this aspect is crucial to ex- 
plain some of the main observations of the experiments 
0, Q and to obtain a complete understanding of BEC 
superradiance. While spatial effects have also been con- 
sidered in Rcfs. 0, Q, the majority of the theoretical 
literature on BEC superradiance is based on spatially- 
independent models where propagation effects are not 
taken into account (see, e.g., Refs. [HlllEEl El) Nev- 
ertheless, the spatially independent description is still 
useful to investigate some basic aspects of the physics un- 
derlying superradiant scattering from BECs, and in this 
paper we will apply both types of models in our analyses. 

The purpose of this paper is to discuss some topics 
relevant to sequential superradiance that were not cov- 
ered in our previous work 0, 0. After briefly review- 
ing the spatially dependent and independent semiclassi- 
cal descriptions of BEC superradiance in Sec. II, we turn 
to the study of sequential scattering in the Bragg regime 
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in Sec. III. It has already been established in Refs. 0,0 
that the characteristic fan-shaped side-mode distribution 
patterns observed in this regime pj can be reproduced 
with the help of the semiclassical models. In this paper 
we would like to discuss several simple, yet important 
physical mechanisms that we think essential for a concep- 
tual understanding of the system behavior in the Bragg 
regime. In particular, they will be useful in explaining 
the emergence of the characteristic fan patterns. More 
specifically, the principles that we will discuss are the fol- 
lowing: first, the concept of "self-adapting" Bragg scat- 
tering, second, the role of the detuning barriers related to 
the transitions between different atomic side modes, and 
third, a directionality in the population transfer between 
side modes that favors the transition along the direction 
of the pump pulse rather then against it. Self-adapting 
Bragg scattering refers to the fact that the optical end- 
fire modes do not contain only a single frequency reso- 
nant with the transition between the BEC at rest and 
the first-order forward side modes. Rather, each tran- 
sition between side modes contributes its own resonant 
frequency component to the endure mode fields. The 
creation of these components ensures the continuation of 
the supcrradiant scattering process. 

A very interesting question regarding the nature of 
the Kapitza-Dirac regime in superradiant scattering was 
raised in Ref. ^lj. Whereas Ref. 0] claimed that 
Kapitza-Dirac scattering is essentially connected to the 
duration of the pump pulse and only occurs at short 
enough times, Ref. [llj stated that the pulse strength 
should be the distinguishing feature. We would like to 
contribute to this discussion by numerically examining 
the regime of strong pulses of longer duration. Our stud- 
ies that are discussed in Sec. IV indicate that for suffi- 
ciently long times, although the system shows a complex 
behavior, one can indeed observe an overall tendency for 
a transition towards the Bragg scattering regime. Apart 
from spontaneous Rayleigh scattering, BEC supcrradi- 
ance can also be initiated by transferring a small amount 
of population into one of the side modes by means of 
Bragg diffraction. This seeding mechanism has been used 
in studies of matter- wave amplification [lH fl3L llij . In 
Sec. V, we investigate how the seeding of a side mode 
modifies the sequential supcrradiant scattering. The pa- 
per concludes with a short summary in Sec. VI. 



II. SEMICLASSICAL THEORY OF 
SUPERRADIANT SCATTERING 



radiant light scattering into the optical endfire modes 
£±(z, t)e y (e~ l ^ tTkz ^ + c.c). As a consequence, the re- 
coiling atoms have well-defined momenta and appear in 
distinct atomic side modes. In the side mode indexed 
(n, m), atoms possess momentum h(nkie x + mke z ), while 
their kinetic energy is given by hu> n , rn w h(n 2 + m 2 )uj r 
with uj r = hkf /2M the atomic recoil frequency. Since the 
pump pulse is far-off resonant, we can adiabatically elim- 
inate the excited atomic state and only need to consider 
atoms in their internal ground state. 

To derive equations of motion for the optical and 
matter- wave fields E(x, i) and ip(x,t), we decompose 
them according to 

^(x,t) = J2 ^ nm ^i ,t K - i ^^ t - nk ' x - mkz \ (i) 

i n v A 

E (+) (x,i) = E^ye-^-^ /2 + £ + (z,t)e y e~ l(uJt - kz) 
+£-{z,t)e y e- l ^ t+kz \ (2) 

with E* + ) denoting the positive-frequency component of 
the electric field, and ip nm (z,t) and £±(z,t) the slowly- 
varying field envelopes along the z axis. The average 
condensate cross section is denoted A, and uj = kc. The 
sum of Eq. JQ| contains only terms with even n + m since 
we only consider atoms in their internal ground state. 
In order to introduce a concise notation, we rescale the 
optical and matter waves as 
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and define the dimcnsionless time t = 2u> r t and length 
£ = kiz. Starting from the full Maxwell-Schrodingcr 
equations for the coupled electric and matter-wave fields 
[la Il6l Il7j , we can use the expansions (JTJ and to- 
gether with the slowly-varying envelope approximation 
to derive the equations of motion for the_enyelope func- 
tions ipnrn and £± . In this way, we find 
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where the coupling constant between matter-wave and 
electric fields is given by 



In this section we briefly summarize our theory of su- 
perradiant Rayleigh scattering developed in Refs. 0,H|. 
As outlined in the Introduction, we consider an elon- 
gated condensate of length L and oriented along the 
z axis that is exposed to a far-detuned pump laser 
pulse £ i (t)e a (e i ( feia; -" !t ) + c.c.)/2 with u x = ck u travel- 
ing in the x direction. Collective effects in the inter- 
action between the atoms and the pulse lead to super- 
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Here, d denotes the atomic dipole moment between 
ground and excited state, 6 the laser detuning from reso- 
nance, and £q the vacuum permittivity. Due to the small 
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size of the condensate, we can neglect retardation effects 
in the propagation of the electric fields. Therefore, we 
can directly express the endfire mode envelopes in terms 
of the atomic side modes, i.e., 



e+(£,r) 



e-(£,T) 



X J — OO i \ 

(n,m) 

X^n m (C,r)tpX+i, m -i^\r), (4b) 
_il r e t{n+rn)T 

xVw(£',T)v>; +1 . m+1 (£',T) (4 C ) 



with x = 2~- From these equations one may clearly see 
how the build-up of the endfire-mode fields at (£, t) is 
driven by the coherences ^nmipn+i m±i- 

The physical contents of Eq. I|4a|) can be interpreted as 
follows. The first term on the right-hand side describes 
the quantum-mechanical dispersion of the envelope func- 
tion, while the second one leads to a spatial translation 
with velocity v m = mftk/M. The remaining terms de- 
scribe photon exchange between one of the endfire modes 
and the laser beam. In particular, through stimulated 
scattering, an atom in a side mode (n, m) can absorb a 
laser photon and deposit it into one of the endfire modes. 
The accompanying recoil transfers the atom into one of 
the side modes (n+1, mil). Alternatively, the atom may 
absorb an endfire-mode photon and emit it into the laser 
beam, thereby ending up in the side mode (n— l,m±' 
This latter process is causes atoms to recoil backward 

In general, superradiant scattering is initiated by 
quantum- mechanical noise, i.e., spontaneous Rayleigh 
scattering from individual condensate atoms 0. Adapt- 
ing the discussion of Refs. |3 DU US regarding "con- 
ventional" superradiance, we can take the effects of 
the initial quantum-mechanical fluctuations into account 
by solving the semiclassical equations of motion with 
stochastic initial conditions (seeds) for the side modes 
0i,±i(£,' r = 0). Since the noise is practically relevant 
for the initial population of those modes only, we can set 
ipnm{ti,'r = 0) = for all other side modes. The con- 
densate -000 is chosen to be in the macroscopic ground 
state at t = 0. In our simulations, we follow the exper- 
imental data of Refs. 0, Q and consider a 87 Rb BEC 
of N = 2 x I0 6 atoms with length L = 200^m and 
cross-section diameter 15/im. The condensate is in the 
Thomas-Fermi regime, s o tha t we can model its wave 
function as ipoo{ z ) = \/n{z) with n(z) = CHL/2) 2 — 
z 2 ]&{L/2 - \z\), C = 3N/4L 3 . As in Refs. 0, ||, the 
first-order atomic side modes are seeded according to 
"0i.ii ( z ?0) = ipoo(z) / "J~N , which corresponds to one de- 
localized atom in each of the modes. Our results are not 
qualitatively influenced by the choice for the shape of 
the seed function. The applied laser pulse is modeled as 
rectangular lasting from t = up to t = tj . 

As mentioned in the Introduction, while the spatially 
dependent description derived above is essential for a 



complete understanding of superradiant scattering, basic 
aspects can already be studied within a simpler spatially 
independent model. The equations of motion for this 
model can formally be derived from Eqs. by dropping 
any spatial dependence and setting ipnm —> C nm / \fk[L 
and e± — > b±/y/kiL. The amplitudes C nm (r) and b±(r) 
can be interpreted as amplitudes of the matter- wave and 
optical fields, respectively, and they obey the equations 
of motion 
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with 7 = 2c/ L. 



III. SEQUENTIAL SCATTERING IN THE 
BRAGG REGIME 



The Kapitza-Dirac and the Bragg regimes of super- 
radiant scattering are characterized by the emergence of 
disctinctive X-shaped and fan-like side-mode distribution 
patterns, respectively. While the appearance of the X 
patterns can concisely be explained in terms of the spa- 
tial localization of the endfire and side modes at the BEC 
edges 0, 13 , the fan patterns show a significantly more 
complex behavior which appears to be determined by an 
interplay of different effects. In this section we would like 
to discuss three concepts which we think to be important 
in understanding the physics underlying sequential scat- 
tering in the Bragg regime. 

(i) Self- adapting Bragg scattering. Following the onset 
of the superradiant process, the optical endfire modes 
contain only one frequency component u> — u>i — 2lj t 
(here we are considering times larger than the recoil time 
so that frequencies can be well defined on the scale of 
Lu r ). However, as more and more side modes become oc- 
cupied, we see from Eqs. (14 b|) and l|4*c|) that additional 
frequency components are created. More specifically, we 
find that the simultaneous presence of population in the 
side modes (n, m) and (n + 1, mil) produces a frequency 
component at u>i — 2(n ± m + 1) u> r . This frequency is res- 
onant with the transition between the two side modes. 

However, to initiate the population transfer to an 
empty higher-order side mode, it is not necessary that 
resonant light is present in the endfire modes. Rather, 
small amounts of population can also be transferred in an 
off-resonant Bragg scattering process involving the pump 
laser, an occupied side mode, and detuned light in the 
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FIG. 1: Frequency spectra |£+(w)[ 2 of emitted endfire mode 
radiation. Calculations are for g — 3.5 x 10 5 s" 1 and pulse 
durations 210 ^s (a), 350 fjs (b), 1100 /is (c); g = 6.4 x 10 5 s" 1 
and duration 750 /xs (d). The insets show gray-scale represen- 
tations of the side mode distributions (as explained in Fig. 01 
and Sec. V) at the end of the pulse. 



endfire mode 1 . As discussed above, as soon as the side 
mode becomes slightly occupied, subsequent stimulated 
Rayleigh scattering into this side mode occurs under the 
production of resonant light since in this way energy con- 
servation is ascertained, and a resonant Bragg scatter- 
ing process ensues. In other words, supcrradiant Bragg 
scattering adapts itself to the occupation of higher-order 
side modes by adding resonant components to the fre- 
quency spectrum of the endfire mode radiation. The self- 
adaptation process is crucial, because otherwise higher- 
order side modes could never be occupied because of the 
increasing detuning barriers. 

In Fig.^ wc illustrate this self-adaptation by calculat- 
ing some spectra |f+(w)| 2 for superradiant light emission 
with different applied laser pulses. Thereby, £{to) denotes 
the temporal Fourier transform of £+(z, t)e~^ UJtTkz \ 
The field envelope £+(z,t) is calculated with the help 
of Eqs. (@J for a rectangular applied laser pulse, and the 
Fourier transform is taken at some fixed zq outside the 
condensate. Figure ^a) shows a situation where only 
the two first-order side modes have been occupied. Wc 
obtain a single broad peak centered at loi — 2io r . As 
soon as the side mode (2,0) becomes populated as well, 
the single peak splits [Fig. E ETb)] . This effect has already 
been analyzed in Refs. [HJ, |2y. The situation of Fig. 
^c) corresponds to a case where the higher-order side 
modes (2, ±2) and (3, ±1) have been populated appre- 



1 Alternatively, the transfer might also be induced by sponta- 
neous Rayleigh scattering from the side-mode atoms, similarly 
to the initial onset of superradiance where the two first-order 
side modes (1, ±1) become seeded from the BEC. It would be of 
interest to compare the relative significance of the two alterna- 
tives. 
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FIG. 2: Grouping of momentum side modes. The transition 
between two side modes is resonant with endfire-mode radi- 
ation of frequency u>i + xu r where x is given by the number 
in the box. The resonance conditions imply detuning barriers 
that have to be overcome in the course of sequential super- 
radiant scattering and lead to a separation of the side modes 
into several groups as indicated by the symbols. 



ciably. We see that another frequency component at 
loi — 6o-v arises. This shows that the endfire modes now 
also contain light resonant with the transition to these 
side modes. In case even higher-order modes become 
populated, two further pronounced peaks appear in the 
spectrum [Fig.^d)]. One is located at &i — 10(J r as ex- 
pected, whereas the other one is shifted slightly from the 
predicted value u>i — lAuj r to toi — 15uv However, we also 
find a stronger background, and the spectrum altogether 
becomes more complex. 

(ii) Detuning barriers. Resonant Bragg scattering be- 
tween the side modes (n, m) and (n + 1, m ± 1) requires 
endfirc-mode radiation of frequency io\ — 2(n ± m + l)uj r . 
Figure displays a selected set of side modes together 
with the resonance frequencies for transitions between 
adjacent modes. The figure shows that the resonance 
conditions imply the existence of several groups of side 
modes that are separated by detuning barriers. In the di- 
agram, these groups are distinguished by different sym- 
bols. The first group consists of the BEC mode (0, 0) 
(large ellipse), the two first-order side modes (1,±1) and 
the mode (2,0) (circles). Immediately after superradiant 
scattering has started up, the endfire modes only contain 
light of frequency u>i — 2uj r . This light is resonant with 
the transitions within the first group. In other words, 
atoms can pass from the BEC to the side mode (2,0) 
without having to overcome a detuning barrier. However, 
all possible transitions from the first group to the adja- 
cent one that consists of the side modes (2, ±2), (3, ±1), 
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and (4, 0) (circles) arc resonant with light of frequency 
lui — 6uj r . This means that the two groups are sepa- 
rated by a detuning barrier. In our discussion of self- 
adapting Bragg scattering, we have explained how the 
initially empty side modes of the second group become 
populated and resonant radiation at frequency u>i — 6uj r 
is produced. However, this process requires some time so 
that the population of the second group sets in only after 
some delay. 

The same reasoning applies to the build-up of popu- 
lation in the subsequent groups of side modes for which 
scattered light at frequency u>i — lOttV (diamond symbols), 
lui — 14u> r (not shown), etc. is required. The height of the 
detuning barrier between adjacent groups stays constant 
at 4cj r . 

(Hi) Directionality of population transfer. The ar- 
guments given in (ii) do not explain why the side- 
mode group shown with hexagonal symbols does not be- 
come populated appreciably in typical experiments in the 
Bragg regime. All the side modes in this group would re- 
quire resonant radiation of frequency u>i + 2uj r to become 
populated from adjacent side modes. This behavior is 
due to a directionality inherent in the population transfer 
between the side modes. This effect is most easily under- 
stood from the spatially- independent model of Eqs. JHJ. 
There, one finds that the population P nm (t) = \C nm (t)\ 2 
of a side mode evolves according to 

dPnm jt) p p p p 

^ r n-l,mil r nm ti+l,m±l r ntn 

+ additional terms (6) 

To derive this relation, we have inserted Eqs. (|5b|) 
and (|5c|) into Eq. JSaJ. The first two terms on the 
right-hand side of (JSJ) arise from the interaction of 
the side mode (n, m) with the resonant components 
e o„ m O n+1 m ±i ana e On-i,m^iL'„ m m 

the endfire mode radiation. They clearly indicate that 
population transfer from the side mode (n, m) preferen- 
tially occurs to the modes (n + l,m ± 1) rather than 
to (n — 1, m ± 1). Physically, this behavior is explained 
by the fact that a forward transfer to the side modes 
(n + 1, m ± 1) amplifies the resonant endfire mode radi- 
ation and in this way the rate of subsequent Bragg scat- 
tering. Conversely, a backward transfer attenuates the 
endfire mode field. The additional terms in Eq. © that 
are not written out are either explicitly oscillating or at 
least do not possess a well-defined sign. They thus do 
not imply a preferred direction in the population trans- 
fer. Nevertheless, they play an important role in the 
overall process since, for example, they can initiate the 
population of empty side modes. 

The concepts explained above can account for the ap- 
pearance of the fan-shaped momentum patterns: The 
atomic side modes can be grouped into sets according 
to the detuning barriers implied by the resonant tran- 
sition frequencies (Fig. |2J). In order to populate a new 
group a detuning barrier has to be overcome which only 
happens after some delay since resonant radiation has to 



be generated. The population of new groups thus oc- 
curs discontinuously in a stepwise process. However, due 
to the directionality of population transfer, only certain 
side-mode groups become occupied (filled symbols in Fig. 
|5J ■ The appearance of the fan patterns then becomes im- 
mediately obvious from the arrangement of the groups in 
Fig. El 

To illustrate the operation of the principles discussed 
above, it is very instructive to study the example pre- 
sented in Figs. 13(a) and (b) of Ref. ||, although the 
full fan pattern does not show up there. These simula- 
tions demonstrate the generic behavior of the side-mode 
dynamics for very weak pulses within the spatially inde- 
pendent model (see also |22j, in particular Fig. 5). Three 
groups of side modes are involved: the first one where 
all modes are temporarily populated, and the two sub- 
sequent ones where only the modes (3, ±1), (4,0) and 
(5, ±1), (6,0), respectively, become occupied. In par- 
ticular, we see that population transfer within a group, 
e.g., from (3,±1) to (4,0), occurs rather quickly once 
the modes (3, ±1) become populated; however, the trans- 
port between groups, e.g., from (4, 0) to (5, ±1) only hap- 
pens after a delay during which resonant radiation is 
built up. For the same reason, population of, e.g., the 
modes (4, ±2) that could be reached from (3,±1) does 
not occur. This transition is detuned so that the atoms 
rather undergo the resonant and fast transfer to (4,0). 
We also point out the temporal regularity in the transi- 
tion process which also persists at longer times that are 
not displayed. For this regularity to appear it is neces- 
sary that the detuning barrier remains the same between 
different side mode groups. Furthermore, the example 
also clearly illustrates the directionality in the transitions 
since the atoms only undergo transition to higher-order 
side modes. In particular, the atoms rather "wait" un- 
til they can overcome a detuning barrier instead of go- 
ing backwards on a resonant transition. Finally, we note 
that full fan patterns are realized for somewhat stronger 
pump pulses. If we consider, e.g., the transition between 
the first and second side-mode group, radiation of fre- 
quency lui — 6uj r can build up before all population has 
moved from the side modes (1, ±1) to (2, 0). In this way, 
the modes (2, ±2) become populated from the (1,±1) 
modes, and the full fan pattern starts to emerge. 

We think that the mechanisms discussed in this sec- 
tion are essential for explaining and understanding se- 
quential scattering in the Bragg regime. Nevertheless, 
the actual behavior found in the numerical simulations is 
rather complex, and the principles discussed here cannot 
account for all observable phenomena. For example, al- 
though we have shown in Ref. |5j that spatial effects are 
not essential for bringing about the fan-shaped distri- 
bution patterns, they may become significant in certain 
situations. In the case of very weak pulses that was dis- 
cussed in the previous paragraph they lead to substantial 
deviations from the behavior found in the spatially inde- 
pendent model Furthermore, they also complicate 
the simple picture conveyed in Eq. (jSJ without, however, 
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FIG. 3: Ratio between backward and forward scattered 
atoms as a function of the pulse duration for g — 5 x 10 6 s _1 
(full curve), 3xl0 6 s _1 (dashed), 10 6 s _1 (dotted), 3 x 10 5 s _1 
(dash-dotted). The inset shows this ratio for fixed field 
strength of the endfire modes e+ = e_ = 4 x 10~ 3 s _1 (full 
curve), and lO - ' 



(dashed) with g = 3 x 10 6 i 



affecting the essential conclusion. The simulations also 
show that population can sometimes move backwards or 
oscillate between side modes. 



IV. TRANSITION BETWEEN KAPITZA-DIRAC 
AND BRAGG REGIME 

The Kapitza-Dirac regime is characterized by the fact 
that individual atom-photon scattering events are not 
rigidly constrained by the condition of energy conserva- 
tion. This means that the atoms may equally well scatter 
laser and endfire-mode photons leading to the simultane- 
ous appearance of backward and forward recoiling atoms. 
In the Bragg regime, the atoms only possess forward- 
directed momenta, i.e., the predominant process is the 
scattering of laser radiation into the endfire modes in ac- 
cordance with the rules of energy conservation. In Ref. 

, an interesting question was raised regarding the dis- 
tinction between the two cases. It was suggested that the 
regime of scattering is determined by the strength of the 
atom-photon coupling g rather than the duration of the 
pump laser as claimed in Q. 

In this section, we would like to further explore this 
issue by means of numerical simulations with strong ap- 
plied pulses of longer duration. Here, "strong" refers 
to a large superradiant gain compared to the recoil fre- 
quency, whereas the scale for measuring the pulse du- 
ration is set by the recoil time. The main results 
of our studies are summarized in Fig. |3J There we 
show the ratio between backward and forward scattered 
atoms, i.e., £ n <o,m P n ,m/ E„>o,m P ™,™ with P «,™ = 
J dz\ip n m (z,t)\ , as a function of the pulse duration for 
different values of the coupling strength. By definition, 
this ratio should be close to I for the system in the 



Kapitza-Dirac regime. Based on these results we can 
draw the following conclusions within our model. 

(i) The Kapitza-Dirac regime is only realized for strong 
and short pulses. With increasing pulse duration, the sys- 
tem enters the Bragg regime which can be seen from the 
sharp drop in the backward-to-forward ratio. We remark 
that in our simulations the X pattern builds up while 
the system is still clearly in the Kapitza-Dirac regime, in 
agreement with experimental results • It is interesting 
to compare our observations regarding the transition be- 
tween Kapitza-Dirac and Bragg regime to a case where 
the amplitudes e±(z,t) are kept constant. This corre- 
sponds to externally apply probe lasers along the BEC 
axis 2 . In this case we observe a behavior displayed in the 
inset of Fig. |3| We see that for sufficiently strong laser 
pulses the ratio between backward and forward scattered 
atoms remains much higher than in the original simu- 
lations and can even exceed one. This shows that the 
dynamic nature of the endfire modes plays a crucial role 
in the transition between the two regimes of scattering. 

(ii) The time that the system spends in the Kapitza- 
Dirac regime is decreased for growing laser intensity. This 
behavior might perhaps qualitatively be explained by the 
obervation that for growing coupling strength g, the sys- 
tem more quickly populates higher-order side modes that 
cover an increasing range A of dctunings. The time scale 
t rsj I/A over which nonrcsonant effects become notice- 
able is therefore shortened. 

(iii) For weak pulses the system is never in the Kapitza- 
Dirac regime, i.e., forward scattering always dominates 
over backward scattering, even at short times. This is in 
agreement with the analytical estimates of Ref. || . 

We have confirmed the conclusions outlined above with 
a complementary analysis of the system dynamics. To 
this end, we have monitored the time evolution of "effec- 
tive rates" for photon emission and absorption to and 
from the optical endfire modes. This provides an al- 
ternative characterization of the two regimes of scatter- 
ing, which of course is not completely independent of 
the one used above. The rates were calculated based 
the population transfer between the atomic side modes 
[the rates of population exchange between adjacent side 
modes can be obtained from Eq. (|4a|l . and, e.g., a net 
transfer from (n, m) to (n + 1, m ± 1) contributes to the 
total photon emission rate]. These studies show that in 
the strong-pulse regime, emission and absorption rates 
are almost equal at short times (Kapitza-Dirac regime). 
Simultaneous with the drop of the curves in Fig. [3 how- 
ever, emission starts to greatly exceed absorption. This 
clearly indicates the transition to the Bragg regime. In- 
terestingly, with increasing pulse duration both emission 
and absorption rates go through a maximum, and sub- 



2 For simplicity, we neglect probe depletion and photon scattering 
between the probe lasers. The intensities of the probe lasers 
are chosen such that the momentum spreads Ap x and Ap z are 
similar to the original simulations. 
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FIG. 4: Atomic side-mode distributions obtained within the 
spatially dependent model and unbalanced seeding. The gray 
level of each square represents the relative probability p nm as 
defined in the text, (a-b) Strong-pulse regime for g — 4.5 x 
10 e s _1 and pulse durations t p = 2.1 (is (a); t p = 3.4/is (b). 
(c-f) Weak-pulse regime for g = 4.5x 10 5 s" 1 and t p = 82.35 /is 
(c); t p = 103.9ms (d); t p = 135.5 ^s (e); t p = 192.6 ^s (f). 



scqucntly the behavior becomes somewhat less clear-cut, 
although emission seems to always be stronger than ab- 
sorption. 



V. SEQUENTIAL SCATTERING WITH 
SEEDED SIDE MODES 



So far we have concentrated on superradiant scatter- 
ing which is initiated by quantum mechanical noise, i.e., 
spontaneous Rayleigh scattering which occurs in the ab- 
sence of any side-mode populations. Experimentally, this 
spontaneous startup of the phenomenon can be strongly 
suppressed by seeding, for instance, one of the first-order 
atomic side modes. Such a preparation can be achieved 
by initially exposing the BEC to a pulsed optical stand- 
ing wave which transfers a small fraction of the atoms 
from the BEC to the desired side- mode [let us say (1,1)] 



by Bragg diffraction |23|. In this case, subsequent appli- 
cation of a pump pulse leads to a very rapid growth of the 
initially occupied side mode due to bosonic enhancement 
whereas the population transfer to the other first-order 
side mode (1,-1), is much slower since it is still initiated 
by quantum mechanical noise. As a result, at the end of 
the pump pulse the main part of the population has been 
transferred from the BEC (0,0) to the side mode (1,1), 
provided that higher-order side modes arc not populated 
H^HEUJI We see therefore that in the case of seeding, 
superradiant scattering essentially leads to amplification 
of the initially prepared weak matter-wave side mode. 
Most importantly, it has been shown that this matter- 
wave amplification is phase coherent [13 , HIH Il4j . 

In this section, we would like to discuss superradiant 
Rayleigh scattering initiated by seeding one of the first- 
order side modes. In particular, our purpose is to explore 
its effect on the distribution of the atomic side modes in 
sequential superradiant scattering. To this end we have 
solved the spatially-dependent Eqs. (4), with various un- 
balanced seed functions ipi.-i(z,0) 7^ ipi,i(z, 0). In Figs. 
01 we present the resulting side-mode distribution pat- 
terns for both strong and weak-pulse regimes with seed 
functions V>i._i(z,0) = ipo t o{z, 0)/V~N and -0i,i( z i 0) = 
V / 50Ch/'o,o (z,0)/ VN- The gray level of each square repre- 



sents the "relative" probability p nm 

1 2 „„A D (max' 



Pn 



j(raax) 



nm with 



IP 

Pnm = J dz\lpnm{z, t)\ 2 and Pnm* — mHj n m )j nm 

In general, the stimulated nature of the startup is ev- 
idenced in two respects. First, the time scales are much 
shorter than in the case of spontaneous startup (e.g., 
compare to Figs. 3 and 11 of Ref. Q). Second, the side 
modes coupled to the initially populated (1,1) side mode 
become populated much faster than other side modes. 
Hence, as depicted in Figs. IHa,b), m t ne strong-pulse 
regime only diagonal side modes (n, n) become popu- 
lated. Such a distribution might be useful for the con- 
struction of a bidirectional correlated atom laser. More 
precisely, it is obvious that if the strong pump pulse is 
followed by a free evolution for a time if , one obtains a se- 
quence of atomic clouds (pulses) propagating in opposite 
directions and separated by distances ~ s/2Hkit{/M, A 
"seeded" realization of the atom laser might have several 
advantages compared to the spontaneous start-up. First, 
one has a much greater control over the time evolution of 
the process since it is not initiated by fluctuations. Sec- 
ond, the effect of potentially detrimental collisions be- 
tween side modes (n, n) and (n, — n) emanating in the 
forward direction (see Fig. 2(a) of Q) is significantly re- 
duced due to the imbalance in the populations. 

In the weak-pulse regime the atomic side-mode distri- 
butions are more involved, in general [see Figs. IHc)- 
(f)]. For relatively short times only side modes in the 
forward upper diagonal (n, n) are populated. This is be- 
cause the field of the 8+ endfire mode initially is too weak 
to induce appreciable population of the (2,0) or other 
off-diagonal side modes (compare with discussion in Sec. 
V.A of 0). Similar to the strong-pulse regime [Figs. 
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Etc)-(f)], we thus typically have new types of patterns in 
the sense that they are clearly asymmetric with respect 
to the laser direction contrary to the case of spontaneous 
startup where they were almost symmetric. The unbal- 
anced seeding has broken this symmetry. However, for 
sufficiently large times bosonic stimulation starts work- 
ing efficiently also for the spontaneously initiated side 
mode (1,-1) and the side modes coupled to it become 
populated. As a result we observe a redistribution of the 
population and the appearance of more symmetric pat- 
terns resembling the well known fan-shape observed in 
the case of spontaneous startup (compare to Figs. 11 of 
Ref. 0). The re-emergence of more symmetric patterns 
appears to be more pronounced than in the strong-pulse 
regime due to the longer time scales involved. 

VI. SUMMARY 

In this paper we have discussed several aspects of 
sequential superradiant scattering from atomic Bosc- 
Einstein condensates. Our treatment was based on a 
scmiclassical description of the process in terms of the 
coupled Maxwcll-Schrodingcr equations for the matter- 
wave and optical fields, developed in Refs. 0,0- First 
of all, we investigated sequential scattering in the Bragg 
regime. We identified three mechanisms which we think 
essential for understanding the system behavior in this 
regime. A first element are the effective detuning barriers 
between the atomic side modes. Although at first sight 
the arrangement of these barriers already seems to im- 
ply the formation of fan-shaped side-mode patterns (see 
Fig-EJ), the actual appearance of these patterns can only 
be explained conclusively by taking two further mech- 
anisms into account: on the one hand, "self-adaptive" 
Bragg scattering by which the system can adapt the fre- 
quency spectrum of the scattered endfire mode radiation 
such that the transitions between atomic side modes re- 
main resonant. If it weren't for this effect, the atoms 
would not be able to occupy higher-order side modes, in 
particular for weaker pump pulses. On the other hand, 
the population transfer between side modes preferentially 
takes place along the direction of the pump laser, rather 
than against it. This aspect makes sure that only certain 
atomic side modes can become populated in the Bragg 



regime. We also pointed out that the system dynamics 
may become quite complex at times and further factors, 
e.g., spatial effects, may have to be taken into consider- 
ation to fully account for the observed behavior. How- 
ever, rather than trying to theoretically analyze all these 
effects in detail, it would be more desirable to have fur- 
ther experimental data available to refine the models and 
to gain more insight into the system. For example, if a 
stepwise population of the different side-mode groups as 
described in Sec. Ill could be verified experimentally, this 
would somewhat favor an interpretation of the superra- 
diant process in terms of atomic diffraction from opti- 
cal gratings, rather than photon diffraction from matter- 
wave gratings (cf. discussion in Q). 

As a second topic, we examined some aspects of the re- 
lationship between the Kapitza-Dirac and Bragg regimes 
of superradiant scattering. Within our model, we find 
that for weak applied pulses, the system is always in the 
Bragg regime, even at times short compared to the recoil 
time. This is of course already expected from analytical 
models of the short-time dynamics In other cases, 

however, analytical descriptions are not feasible, and nu- 
merical methods have to be applied. In particular, when 
investigating stronger pulses, wc find that the system al- 
ways makes a transition from the Kapitza-Dirac to the 
Bragg regime for increasing pulse duration. Interestingly, 
the time spent in the Kapitza-Dirac regime seems to de- 
crease for growing pulse duration. 

Finally, we also investigated how the initial seeding of 
one of the atomic side modes by way of Bragg diffraction 
influences the sequential scattering process. Questions of 
this kind can be expected to be of interest in connection 
with studies of matter- wave amplification 0, fl3L Il4j . 
We presented typical examples for atomic side mode dis- 
tributions in the weak- and strong-pulse regimes. In 
particular, we find new patterns in the sense that they 
are clearly asymmetric with respect to the laser direc- 
tion contrary to the case of spontaneous startup. This 
asymmetry is of course reflective of the preparation of 
the initial state. The extent of asymmetry however, can 
vary appreciably in the course of time, especially in the 
weak-pulse regime. 

GMN acknowledges support by "Pythagoras II" of the 
EPEAEK research programme. 
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